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ABSTRACT 


An optical digital communication system requires that an accurate 
clock signal be available at the receiver for proper synchronization 
with the ' transmitted signal. Phase synchronization is especially 
critical in M-ary pulse position modulation (PPM) systems where the 
optimum decision scheme is an energy detector which compares the energy 
in each of M time slots to decide which of M possible words was sent. 
Timing errors cause energy spillover into adjacent time slots (a form of 
intersymbol interference) so that only a portion of the signal energy 
may be attributed to the correct time slot. This effect decreases the 
effective signal, increases the effective noise, and increases the 
probability of error. 

In this report, we simulate a timing subsystem for a 
satellite-to-satellite optical PPM communication link. The receiver 
employs direct photodetection, preprocessing of the detected signal, and 
a phase-locked loop for timing synchronization. The variance of the 
relative phase error is examined under varying signal strength 
conditions as an indication of loop performance, and simulation results 
are compared to theoretical calculations. 
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CHAPTER 1 
INTRODUCTION 


In an optical digital communication system, the laser is modulated 
using an appropriate signaling format. In pulse position modulation 
(PPM) systems, an accurate clock signal must be available at the 
receiver to provide proper synchronization with the transmitted signal. 
In some systems, timing information is sent to the receiver in the form 
of a special timing signal which is superimposed on the transmitted data 
waveform. This timing signal is identified and separated at the 
receiver and used to ensure time synchronization. Alternatively, a 
receiver could contain a subsystem to extract timing directly from the 
data stream and no special timing waveform would be required, resulting 
in a simpler and more power efficient transmitter. A synchronous 
subsystem usually contains a local oscillator tuned near the expected 
frequency of the incoming signal. Correct frequency tuning of the local 
oscillator is essential for synchronous reception. However, the phase 
of the local oscillator must also agree with that of the incoming 
signal. Phase synchronization is especially critical in M-ary PPM 
systems where the optimum decision scheme is an energy detector which 
compares the energy in each of M time slots in a frame to decide which 
of M possible words was sent. A timing error causes energy spillover 
into adjacent time slots (a form of intersymbol interference) so that 
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only a portion of the signal energy may be attributed to the correct 
time slot. This effect increases the probability of error. The system 
must also have the capability to track frequency displacements due to an 
actual frequency difference between the receiver and transmitter, a 
drift of the receiver oscillator, or Doppler shift, due to motion 
between the transmitter and receiver, as in satellite tracking. 

The effects of timing errors in optical digital systems have been 
extensively investigated [lj. Many methods have been proposed and 
employed for tracking systems including pulse-edge tracking, early-late 
gate detection, and maximum a posteriori (MAP) estimation. Phase-locked 
loops have also been studied for use in both optical on-off keying and 
optical sinusoidal modulation schemes [2], [3], 

In this report, we simulate a timing subsystem for a 
satellite-to-satellite optical PPM communication link. The receiver 
employs direct photodetection, preprocessing of the optical signal, and 
a phase-locked loop for timing synchronization. The photodetector 
output is modeled as a filtered, doubly stochastic Poisson shot noise 
process. The variance of the relative phase error is examined under 
varying signal strength conditions as an indication of loop performance, 
and simulation results are compared to theoretical relations. 
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CHAPTER 2 

DETECTOR" MODEL AND SIGNAL PREPROCESSING 


A diagram of the proposed receiver timing subsystem is given in 
Figure 2.1. An optical signal from a free space link or optical fiber 
is incident on a photodetector which produces an electrical current in 
response to the incident photons. The photodetector output is 
preprocessed and applied to the input of a phase-locked loop (PLL) which 
tracks the signal and supplies timing information to the decision 
section of the receiver for proper decoding of the data sequence. 

2 . 1 Pulse Position Modulation 

The signaling method under consideration is M-ary pulse position 
modulation (PPM). In PPM, a single pulse is sent in one of M time slots 
comprising a word or frame. Therefore, M distinct messages can be sent 
during each frame. In the absence of timing errors, the optimum 
detection method for optical PPM is an energy detector or photon 
counter. In this decoding scheme, the receiver compares the number of 
counts in each of the M time slots and selects the word corresponding to 
the maximum count. If M=2 , then the PPM signaling format can be used 
for block encoding, where M=2 represents a binary word of k bits. If 
each time slot is T seconds long, the system transmits data at a rate of 
R = l(log 2 M)/MT) bits/sec. 
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Higher-order PPM systems demonstrate superior performance and 
require less signal power for a given probability of error than other 
digital signaling methods. However, use of PPM carries stringent timing 
requirements. Timing offsets during decoding cause only a portion of 
the signal energy to be included in the signal count for the correct 
time slot. Some energy will spill over into adjacent slots causing 
intersymbol interference. This spillover decreases the effective signal 
power, increases the effective noise, and increases the probability of 
error. 

The performance criterion for block encoded systems is the average 
probability of making a word error, PWE. From [4], the PWE for a block 
encoded PPM signaling format with no timing error is 
exp[-(K + MK,)] 

PWE = 1 jj ° l Pos(k,K s + 1^) 

k=l 


k-1 

l Pos(j.IC) 

j~0 ^ 


M-l 


fl 


Ma 


[(1 + a) M - 1] 


( 2 . 1 ) 


where 


a 



Pos(k,m ) 
v 


A 

S3 


(m ) k 
v 

k! 


exp(-m ) 
v 


and 

K g = average signal count 

K, = average background count 
b 

M = time slots/word. 


( 2 . 2 ) 


The probability of bit error, PBE, is related to the word error 
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probability by 


PBE 


f 

k M 
2 lM- 1 


PWE . 


(2.3) 


Equation (2.3) is useful for comparing the performance of block encoded 
PPM to other digital signaling techniques. 

2.2 Photodetector Model 

The output current of a photodetector, i D (t), consists of a signal 
current component, i (t), due to the incident optical power on the 

» s 



that 


i D (t) - i s (t) + i g (t). (2.4) 

The signal component i g (t) of (2.4) is a non-stationary , inhomogeneous 
shot noise process represented by [4] 

N(t) 

i (t) = l G, w( t - t.) (2.5) 

S J-l 2 J 


where 


J 

N(t) 

G. 

J 

w(t) 


arrival time of the photon 

the total number of arrival photons during (0,t) 
random, statistically independent detector gains 
response of the detector to a single electron. 


In response to each arriving photon, the detector produces a primary 
electron-hole pair. The photon arrival times and thus the electron 
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release times can be modeled as a filtered, doubly stochastic Poisson 
process with arrival rate 

OO 

V° = X o + X s £ p(t - kT w ' C k T W /M) (2 * 6) 

k=-<*> 


where the received optical pulse shape, 
square pulse of duration T w /M shown 
signal counting rate defined by 


A = 
s 


n,P„ 

d o 

hf 


p(t), is a unit-amplitude, 
in Figure 2.2, where A g is the 


(2.7) 


and 

P q = optical power of the received signal 
n d = quantum efficiency of the detector 
h = Planck's constant 

f = the optical frequency of the signal 

M = the number of time slots in a word 

c^ = (0, 1 , 2 , . . . , M} is the random data sequence 

A q = counting rate due to dark current 

T = word width 
w 

T /M = time slot width, 
w 

Note that A^(t) is the average rate of the photon arrival times. In any 
given time period of T seconds, the probability that exactly N counts 
are detected is given by 


p lN, (t o,t 0 + T)] 


N -p 
U e 

N! 


( 2 . 8 ) 


where 
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t +T 


r 0 

v = J 


t 

o 


^(t) dt 


(2.9) 


2.3 Signal Preprocessing 

Correct operation of the PLL requires that a dominating frequency 
component, close to that of the local oscillator in the PLL, be present 
in the input signal. This frequency component is what the PLL "locks" 
onto, enabling it to precisely track the input signal. In analog radio 
communication, for example, the tracked component is the modulating 
waveform frequency. If we expect successful operation of the PLL under 
the proposed operating conditions, we must ensure that the PLL input 
contains an adequate frequency component. 

First consider the output from the photodetector described in 
Section 2.2. Assuming a constant detector gain, a square optical signal 
pulse, p(t), and neglecting thermal noise, the expectation of the 
photodetector output conditioned on the random data sequence, c^, is 

E£i D (t)|c k } - w(t) * ^(t) (2.10) 

where A^t) is the photocount rate given in (2.6) and w(t) is the 
detector impulse response. The unconditioned expectation of the 
detector output is 

E{i D (t)> - w(t) * <A R (t)>' (2.11) 

where 

00 M- 1 . 

< - A + * l l jf p(t - ki„ - jt„/M) . 

k=-00 j=0 


( 2 . 12 ) 
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The second term of (2.12) is the signal pulse shape averaged over the M 
slots in a word. Since the received signal pulse shape, p(t), is 
square, the second term of (2.12) is a constant. Therefore, (2.12) 
reduces to 

X 

<X R (t) > = A o (2.13) 

and 

E {ip( t) } = constant (2.14) 

which implies 

F t E (ijj( t) }] = 2 ir( constant) <S(u>). (2.15) 

From this analysis, we conclude that, for a square pulse shape, there is 

no trackable frequency component in an expected detector output signal, 

0 

only a dc term is present. Modeling p(t) as a square pulse represents a 
worst case situation. If p(t) has a shape other than square, (2.13) 
would no longer hold and other frequency components would appear in the 
output spectrum. However, further analysis would be required to 
determine if the frequency components generated would be of sufficient 
amplitude for the PLL to track. 

We now consider a form of signal processing which will result in a 
trackable frequency component in the spectrum of the signal before it is 
applied to the PLL input. Figure 2.3 is a block diagram of the proposed 
detector and signal preprocessor model. The use of a linear operation 
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will result only in a convolution in the time domain or multiplication 
in the frequency domain which will change nothing but the magnitude of 
the output spectrum. A nonlinear operation is needed to produce a 
frequency component at the slot rate. A square-law device was selected 
to produce the nonlinear effect because of its simplicity and ease of 
implementation and analysis. 

Consider the expected value of the square of the detector output 
current. From Reference [4] we have 

E{ip(t)} = <[h(t) * yt)] 2 > + h 2 (t) * <yt» (2.16) 

where h(t) is the combined impulse response of the detector and an 
additional linear filter contained in the preprocessor. Again, we 
assume a square, unit-amplitude signal pulse, p(t). Since p(t) is a 
square pulse, <yt)> is again a constant and thus the second term of 
(2.16) is also a constant. However, the first term in (2.16) can be 
written as 

00 00 

h(t) * yt) = X o / h(t) dt + X g l g(t - kT w - c k T w /M) (2.17) 

_oo k=-°° 

where 

g(t) = h(t) * p(t). (2.18) 

For computational convenience let 

OO 

/ h(t) dt - 1 . (2.19) 

— 00 


In this case, we obtain 
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o 2 A a 

<[h(t) * yt)j 2 > = a o+ — g-s. 

OO 00 

+ X s l l <g(t " - c^T/M) (2.20) 

j=-oo k=-“ 


*g(t " jT - c.T /M)> . 
w j w 

The first two terms of (2.20) are constants. Assuming the pulses 

occupying the time slots are narrow and no significant overlap occurs, 
the cross terms of the double summation in (2.20) are negligible. The 

more slots per word, the better this approximation will be. Use of this 

approximation implies that the only significant terms in the summation 
are those for which j=k, and we obtain 

A s l l <•> * X 2 s l <g 2 (t - kT w - c k T w /M)> (2.21) 

J = — OO k = “ oo k=-<» 

so that 


** l <g 2 (t - kT w - c k T w /M)> = X l l ¥ g2(t ' nT w /M) • (2 - 22) 

k=-« n=-<*> 

After substituting (2.22), (2.21), and (2.20) into (2.16), we obtain 


A 

00 

o 2 A \ 


OO 

, , s 

/ h 2 (t) dt + 

— OO 

2 so 


/ h(t) dt 

- — OO 

X o + TT 

_ o + M 



+ 




g 2 (t - kT /M) . 

w 


(2.23) 
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Since all but the last term are constants, (2.23) is equivalently 
expressed as 


E{i£(t)} a K [ h( t ) * p(t)] 2 * l 6(t - kT w /M) (2.24) 

k=-® 


where K denotes the sum of the constant terms of (2.23). 

Taking the Fourier transform of (2.24) yields 

„ ,2 

2 ^ 

[E{ip(t)}] = 2 itK5(w) + [H(oj)P(m)j * [H( co)P( co) J 

w 


(2.25) 


x l 6( <d - k2 irM/T ) 
k=-» w 


where H( u>) and P( w) are the Fourier transforms of h(t) and p(t), 
respectively. P( u>) is simply the transform of the transmitted square 
pulse and is given by 


w 


(JT 


> nr -i« si»<>Jr w /2M) -jjj a 

p <») - / e ■ dc 7z e 

0 w 


(2.26) 


From (2.25) it is apparent that we have spectral frequency components at 


2 xM 

a) = — — n 
n T 

w 


(2.27) 


When n=l, corresponds to the time slot transmission rate. H( w) 
should be chosen to maximize the signal component at and to minimize 
the other background noise components. 

The next step is to put a constraint on the choice of the combined 
receiver filter impulse response. First, the detector output can be 
represented by 


i(t) = i(t) + Ai(t) 


(2.28) 
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where the subscripts are dropped for convenience, and i(t) = E{i(t)}. 
The squared output of the detector is 

i 2 (t) = i ( t ) 2 + 2i (i) Ai(t) + Ai 2 (t) . (2.29) 


Taking expectations on both sides of (2.29) and assuming E{Ai(t)} = 0 
give 

i 2 (t) = TTtT 2 + Ai 2 (t) . (2.30) 

Substituting (2.30) into (2.29) yields 

i 2 (t) =. i 2 (t) + 2i(tT Ai( t) + [ Ai 2 (t) - Ai 2 (t)] . (2.31) 

In comparision to the first two terms of (2.31), the third term, 
2 2 

Ai (t) - Ai (t) , is considered small and can be neglected for a first 
approximation, leaving 

i 2 (t) = i 2 (t) + 2UtT Ai(t) . (2.32) 


To examine the power density spectrum of (2.32), we first calculate 
the autocorrelation function 


<i 2 (t ) ,i 2 (t_)> = i 2 (t ) i 2 (t ) + 4i(t, ) i(t ? ) 

1 2 1 1 (2.33) 

x [i(tj) - i(t 1 )][i(t 2 ) - i ( t 2 ) J 


The first term of (2.33) is the expected output signal from the squaring 
loop which contains a nonzero frequency component; the second term is 
noise. To minimize the noise term, remember 
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i(t) = h(t) * <X R (t)>. 


(2.34) 


So if H( to) is chosen such that 


H(0) = 0 


(2.35) 


the first noise term in (2.33) vanishes. Many functions satisfy (2.35) 
and an optimal choice of H( oi) is beyond the scope of this report. An 
intuitive choice is the bandlimited derivative function of Figure 2.4a, 
with impulse response, h(t), plotted in Figure 2.4b. It can be shown 
that the optimum bandwidth, B q , is equal to the time slot frequency. A 
sinusoidal function also satisfies (2.35). In the simulation model 


h(t) ■ sin 


2 irM 


w 




(2.36) 


which is plotted in Figure 2.5a. The Fourier transform of h(t) is given 
by 


H( oi) = 


irT 

w 

Sln t-H 

( 2 itM-i 

1“ " T J 

w 

T 

sln tzff 

r , 2 ttM i 

1^ + — JJ 

w 

jM 

T 

w r 

2M 

2ttM>, 

T > 
w 

T 

w 

2M 

r , 2tiM>| 

l" 5 + T — J 

w 


(2.37) 


and plotted in Figure 2.5b. This filter was used because of its easier, 
more efficient implementation in software. Preprocessing the detector 
output as described above will generate a spectral component at the time 
slot frequency which can now be tracked by the PLL. 
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Figure 2.5a. Impulse response of sinusoidal filter used in the simulation 
model. 



Figure 2.5b. Frequency response of the sinusoidal filter used in the 

simulation model. This filter approximates the ideal low- 
pass differentiator plotted in Figure 2.4a. 
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CHAPTER 3 

DESIGN AND ANALYSIS OF THE PHASE-LOCKED LOOP 

3.1 Mathematical Model of the Phase-Locked Loop 

A block diagram of a phase-locked loop (PLL) is given in Figure 3.1. 
In the basic operation of the PLL, a phase detector measures the phase 
difference between its two inputs and then outputs a voltage 
proportional to this difference, v d (t). The difference voltage is then 
passed through a loop filter and fed into a voltage controlled 
oscillator (VCO) which changes its frequency in a direction to minimize 
the phase difference between the VCO output and the input signal. The 
loop is said to be "locked” when the control voltage applied to the VCO 
keeps the frequency of the VCO equal to the average frequency of the 
incoming signal. 

For the analysis of the PLL, assume that the input signal has phase 

0 s (t), the VCO output has phase 9 q C t ) , and that the loop is locked. The 

phase detector is modeled as a perfect multiplier of gain K m whose 

output, over a specified operating range, is a direct voltage, v^Ct), 

that is a function of the phase difference 9 (t) - 9 (t). Under locked 

s o 

conditions the phase difference will be small and the phase detector 
output is l 5] 

V s> ' K d le s (t) - 6 o (t)J 


(3.1) 
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where 


K = K A A 
d m o s 


(3.2) 


is the phase detector gain and A q and A g are the amplitudes of the VCO 
output and input signals, respectively. 

To reject high frequency signal components, the phase detector 
voltage is filtered by a low-pass loop filter of transfer function F(s). 
The resulting control voltage is then applied to the VCO. 

The VCO has a free-running frequency and an instantaneous 
frequency shift proportional to the control voltage input, AaFK Q v c (t), 
where K q is the VCO gain constant and v c (t) is the filter output 
voltage. The VCO output frequency is written as 

cOq = o) f + = <o f + K Q v c (t) . (3.3) 


Noting that frequency is the derivative of phase, we see that 


de o (t) 

dt 


= K v (t) 
o c 


(3.4) 


Using the Laplace tranforms of (3.1) and (3.4) we arrive at the 
linearized mathematical PLL model of Figure 3.2. The loop transfer 
function of Figure 3.2 is 


9 (t) K K , F ( s ) 

o o d 

e (t) s + K k.f(s) 

s o d 


(3.5) 


For a further analysis we need to specify F(s). Using a simple 
low-pass filter, F(s) is of the form 
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F(s) 


1 

S T + 1 


(3.6) 


where 1/t is the 3 db cutoff frequency of the filter. Substituting 
(3.6) into (3.5), the loop transfer function becomes 
KK./t 

H(s) =-r — . (3.7) 

s +s/t+Kk,/t 

O a 


From control theory let us define £ as the damping ratio and as the 
natural undamped frequency. The denominator of (3.7) is the 
characteristic equation of the system that determines the system 
response and is defined by 


A = s 2 + 2?(^s + 


which implies 


% 


K K , 
o d 


1/2 


and 


5 = 


tK K . 
o d 


1/2 


(3.8) 


(3.9) 


(3.10) 


The selection of the damping ratio is of particular importance. With 
too large a value of £, the system response is overdamped, resulting in 
a very long response time to an input. With too low a value, the system 
is underdamped, causing oscillation about the steady-state value and 
requiring a long time for a settled final response. 

Two final design parameters to consider are the lock-in range and 
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hold-in range, Au^ and Au^, respectively. When the difference between 
the free-running VCO frequency and the average incoming frequency is 
less than the 3 db loop bandwidth, the loop will lock up almost 
immediately without skipping cycles. The maximum frequency for which 
this is possible is the lock-in range and is roughly equivalent to the 
loop bandwidth [5]. The range of frequency difference over which a loop 
will maintain phase lock is the hold-in range and is defined by [5] 

Aa k = ±K o K d F(s) * (3.11) 

More detailed discussion and derivation of PLL design parameters can be 
found in [5] and [6]. 

3.2 Selection of Design Parameters 

Successful PLL design represents a compromise between performance 
features. A wide loop bandwidth will result in a fast acquisition and 
lock-up time as well as an extended lock-in range, but will admit more 
noise, may have a smaller hold-in range, and be more sensitive to small 
changes in the input signal frequency. Alternatively, a narrow loop 
bandwidth insures good noise performance and will tolerate a larger 
amount of input frequency change without losing lock, but will have 
difficulty initially acquiring phase-lock. Advanced techniques such as 
frequency sweeping and loop bandwidth changing have been developed for 
acquiring lock even with very narrow bandwidths , but they are beyond the 
scope of this report. A thorough treatment may be found in [6], 
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3.2.1 Simulated PLL parameters 

The simulations described in this report used one-second time slots 

which correspond to an expected input frequency (slot rate) of 

2ir rad/sec. The VCO free-running frequency was selected as 

6.18 rad/sec. The difference of 0.103 rad/sec represents a 1.6% 

deviation from the expected input frequency, and is used to account for 

uncertainties about the exact frequency of the received signal. Such 

uncertainties are commonly caused by Doppler shifts, drift in the 

transmitted clock, or drift of the receiver local oscillator. The phase 

detector gain, VCO gain constant, loop filter bandwidth, and damping 

ratio were chosen as: 

K d = 0.3333 V/rad 

K = 1.0 rad/sec/V 

o 

1/t= 0.3333 rad/sec 

£ = 0.5 

which yield the parameters: 

Au^ = 0.3333 rad/sec 

Au^ = 0.42 rad/sec. 

It is important to recall from (2.2) that the phase detector gain 
is dependent on the amplitude of the input signal, A g . In fact, we are 
implicitly assuming that the average magnitude of A g is equal to one. 
If the value of changes due to a change in A g , the design parameters 
of the PLL also change. For example, if A g increases from 1 to 10, the 
3 db loop bandwidth increases from 0.42 to 1.6 rad/sec causing the loop 
to reject less noise than originally intended. Realistically, the input 
signal will vary over a substantial range of values and the square-law 
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device in the preprocessor compounds the problem. To compensate for 
varying signal amplitudes, an attenuator is placed between the squaring 
loop and the PLL input. This device will adjust the preprocessor output 
signal to insure the average PLL input signal will not exceed unity. 

To choose the attenuation constant, we first must find the expected 
value of the amplitude out of the preprocessing components. As in 
Chapter 2, let g(t) denote the convolution of the input pulse shape, 
p(t), and the combined receiver impulse response, h(t). Let each pulse 
have amplitude P q and width T=T w /M sec. From (2.33), the receiver 
impulse response is 


h(t) = sin 




(3.12) 


Convolving h(t) with p(t) yields 


g(t) = 


P T 
o 

2 7T 


P T 
o 

2 ir 


1 2 7T 

1 - cos — t 


i 2 it 

1 - cos — t 


-T < t < 0 


0 < t < T 


(3.13) 


Figure 3.3a shows that as a result of the convolution operation, g(t) 
has duration 2T. After passing g(t) through the squaring loop, we 
obtain 

3P^T^ P^T^ P^T^ 

g 2 (t) = — V- - -V cos 4 1 t +-V cos T 1 t, - T < t<T , (3.14) 

8it 2/ 1 8 it 


which consists of two identical, positive waveforms each of duration T, 

as shown in Figure 3.3b. Thus, the total amplitude of the fundamental 

2 

frequency component of g (t) is 




Figure 3.3a. Result of convolving a sinusoidal receiver impulse response, 
h(t) , with a square signal pulse, p(t). 



2 2 


Figure 3.3b. Result of passing g(t) of Figure 3.3a through the squaring 
loop. 
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= 2 


2 2 
PI 
o 


2 u 


(3.15) 


Averaging (3.15) over the M time slots in a word and substituting 
T = T w /M give us an amplitude coefficient of 


a 


1 


2 2 
P T 


o w 
M 3 ir 2 


(3.16) 


So multiplying the PLL input by the attenuation constant, a, defined by 


1 



mV 

p 2 t 2 

o w 


(3.17) 


limits the average PLL input signal to unity. This completes the design 
and analysis of the receiver timing subsection. 
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CHAPTER 4 

COMMUNICATION SYSTEM LINK 


Of importance in any practical communication link is the amount of 
signal power received. The larger the received signal, the better we 
can determine the original message and the lower the probability of 
making an error. In this analysis, we are ultimately interested in the 
number of photocounts per signal pulse or time slot received by the 
photodetector. Specifically, we now consider optical transmission 
between two satellites, one in a geosynchronous orbit and the other in a 
lower earth orbit. 


4 . 1 Transmitter Design 

The purpose of an optical transmitter is to focus as much of the 
modulated signal power (led or laser) as possible onto the receiving 
telescope. For a laser transmission system, the transmitter usually 
consists of a beam expander which is used to reduce the beam divergence. 
The beam expander may consist of reflecting or refracting optics 
depending on cost, weight, and wavelength of the source. Figure 4.1 
shows two beam expander configurations. 

The output of the beam expander can be considered as a uniformly 
illuminated circular diffracting aperture of diameter d^,. The amplitude 
distribution of the Fraunhofer (far-field) diffraction pattern is [7] 




CASSEGRAINIAN TELESCOPE 


Figure 4.1. Typical bead expander configurations. 
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, k<£[2J 1 (kd T r Hi 
U(r Q ) = exp(jkz) exp(jkr o /2z) 


where 

r Q = radial coordinate in the far-field plane 
d^ = diameter of the transmitting aperture 
k = wave number of the source 
z = transmission distance 


and Jj(x) is a first-order Bessel function. Thus the far-field 
intensity distribution is 


Ur Q ) - U(r Q ) U*(r Q ) = 


2J 1 (kd T r Q /2z) 

kd T r o /2z 


(4.2) 


where I T is the transmitted intensity. This intensity distribution is 
commonly known as the Airy pattern and has well-defined zeros wherever 
the quantity 

= 0 . (4.3) 


2Jj(kdj,r Q /2z) 

kd T r Q /2z 


For large distances, the receiver aperture is small compared to the 
size of the Airy disk (the circular area contained within the first 
minimum of (4.3)). Therefore, we make the approximation for an aperture 
of diameter d^, centered at a point r Q in the receiving plane, that the 
total optical power, P r , intercepted by the aperture is 

,d 2 

P, =-T- * A r « r 0 ) (4 ' 4> 

where A r = ird r /4 is the receiver aperture area. Also, assuming uniform 
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illumination of the transmitting aperture, the transmitter laser power 
is related to the transmitted intensity by 


.2 

P L ~ ~T~ X T = A T I T 


where A^, = ird^/ 4 is the transmitter aperture area. 
(4.4), and (4.5), we write 

} 2 

P = 
r 


Vr 

2J 1 (ird T r Q /zX) 

z 2 X 2 

ird T r o /zX 


(4.5) 

Now using (4.2), 

(4.6) 


where we have used k=2Tt/X, where X is the wavelength of the source. 

We now consider optical power losses in the transmitter. First, no 
optical surface will transmit all of the incident optical power. 
Compounding this problem is the fact that there is generally more than 
one lens or mirror, and the combined effects are multiplicative. 
Second, since the laser cross-section can be mathematically described as 
a Gaussian beam which is infinite in extent, we must account for 
truncation of the expanded beam by the finite diameter transmission 
aperture. Finally, if a Cassegrainian or similar type of reflecting 
telescope is used, the secondary mirror obstructs the beam and reduces 
the overall transmissivity of the telescope. In some cases, the 
secondary mirror mounting supports may also contribute to obscuration of 
the beam. A detailed analysis of beam truncation and obscuration 
effects can be found in [8] and [9]. 

For this analysis, we account for all the transmitter losses by 
defining a constant, n^,, which is the overall efficiency (transmittance) 
of the transmitting telescope. Thus, (4.6) becomes 
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P 

r 


VA 

2,2 

z X 


2 J ^ ( 7rd T r o /2z) 
nd T r o /2z 


(4.7) 


4.2 Receiver Design 

The objective of the optical receiver is to focus the light gathered 
by the aperture onto a photodetector while rejecting as much background 
radiation as possible. A typical receiver is shown in Figure 4.2. The 
field stop iris is used to narrow the field of view of the receiver so 
radiation from areas adjacent to the transmitter is ignored. The 
interference filter eliminates wavelengths which are not the same as 
that of the source. In addition to the optical losses mentioned in 
Section 4.1, the receiver has the additional loss associated with the 
interference filter. Narrow band filters have small transraittances 
whereas wider band filters have higher transmittances. For our 

analysis, as in the transmitter, we define a constant, r^, to be the 
overall receiver efficiency to account for all losses. 

The receiver focuses the optical power onto a photodetector such as 
an avalanche or PIN photodiode. In turn, the incident photons cause the 
generation of electron-hole pairs which flow to produce current in the 
device. The detector quantum efficiency, is defined as 

Number of electron hole pairs generated 
^d Number of incident photons * (-+*8) 

A practical photodiode may have a quantum efficiency of between 30% and 
90%. However, a high quantum efficiency requires a thick depletion 


region of material which results in a slow response time. 


A fast 
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response photodiode, necessary for high data rate communication systems, 
will typically have a much lower quantum efficiency due to the tradeoff 
between response time and efficiency. 

We are interested in determining the number of photocounts per 
signal pulse generated by the photodiode. Letting D denote the data 
rate of the signal, we write 


Dhc r 


(4.9) 


where 

<N> = expected number of photocounts per signal pulse 

h = Planck's constant 

c = speed of light. 

Combining (4.9) and (4.7) gives us 

12 


<N> = 


n d n r n T A T A r 

r 

2J^ ird T r o /zX) 

Dhcz^ X 

w^r /zX 

. 


(4.10) 


Equation (4.10) is the total link equation for the system. Note that 
(4.10) is maximized if r Q = 0 which corresponds to zero pointing error 


between the transmitter and receiver. So 
n d n r n T A T A r „ 


<N> 


max 


Dhcz^ X 


(4.11) 


is the maximum number of photocounts per signal pulse. 


4.3 Background Radiation 


For two satellites in an earth orbit, background radiation (noise) 
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sources include the stars, planets, earth, moon, and sun. Reception 

with the sun directly in the background is usually not attainable 

because of its high brightness, and will not be considered here. A 

useful radiation measurement is the spectral radiant emittance, W(X), 

which is the radiant power at a wavelength X emitted into a hemisphere 

per unit area of source in the hemisphere. Assuming a Lambertian source 

of diameter d at distance z from the receiver, the total amount of 
s 

background radiation power incident on the detector, P is 

D 


P B = 


TTT^ X_.d_d 


V 

i s r 


W(X) 


16z 


(4.12) 


where 


n r = receiver efficiency 

d^ *> source aperture diameter (m) 

d^ = receiver aperture diameter (m) 

z = distance between source and receiver (m) 

X 1 = interference filter bandwidth (A). 

Equation (4.12) assumes that the whole source is contained in the 
receiver field of view. If the solid angle of the receiver field of 
view, Jl r , is less than the solid angle subtended by the source at the 
receiver, , the background radiation power is given to a good 
approximation by [10] 


P B = 


2 2 

*Vi d s d r 

V 

16z 2 



W( X) 


(4.13) 


The solid angle for the receiver (assuming small angles) is 
' ttB 2 

n = 

r 4 


(4.14) 
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where 0 is the planar angle of the receiver field of view. 

(4.14), (4.13) and 



Using 


(4.15) 



(4.16) 


which is the formula for an extended source filling the receiver field 
of view, as is the case when a satellite is pointed towards the earth 
for signal reception. 

For smaller sources such as stars, it is usually more convenient to 
make background radiation measurements in terms of the spectral 
irradiance, H( X) , which is the power per unit wavelength interval 
incident on a unit area of the receiver [10]. Therefore, if a 
background radiation measurement is given in terms of H( X) , then the 


background radiation power gathered by the receiver is 


"Wr 


H(X) . 


(4.17) 


Of primary concern is the background radiation from the earth due to 
reflection of the sun. Figure 4.3 shows the spectral radiant emittance 
of the earth versus wavelength. Note that for the wavelengths of 
interest (between 0.8 and 0.9 iim) , W( X) is near its peak value. As 
indicated in Figure 4.3, this measurement is valid when there is no 
cloud cover in the earth's atmosphere. Radiance from sunlit clouds is 



Spectral radiant emittance, watts per cm 2 , micron 
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approximately one order of magnitude larger and, as we will see in 
calculations to follow, it represents the strongest, worst case 
background noise source. Figure 4.4 shows the spectral irradiance for 
the planets and moon, and Figure 4.5 shows the spectral irradiance for 
the brighter stars immediately outside the earth’s atmosphere. 
Figures 4.3 and 4.4 show that a full moon exibits the largest background 
emission next to that of the earth. Because of the close proximity of 
the satellite to the moon, the use of (4.17) would be incorrect since it 
assumes the whole source is in the field of view, when actually only a 
portion of the moon would be in the satellite's field of view. Assuming 
that the moon appears as a uniformly illuminated Lambertian surface, we 
have the relation [10] 

W( X) = H( X) . ' (4.18) 


Substituting (4.18) and (4.15) into (4.16) yields 


2 2 2 
"h A. 9d z 

P = — - a— - — H( A) 

B 4d 2 

s 


which is the background power gathered by the receiver when the source 
fills the field of view and the spectral irradiance is the given 
measurement parameter. We can convert the background radiation power 
gathered by the receiver aperture into background counts/time slot 
output by the detector using 


<n> b 


n. A 

— — P 
DhcM B 


(4.20) 


where 



10-4 
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Figure 4.4. Calculated planetary and lunar spectral irradiance outside the terrestrial 
atmosphere. (From W. K. Pratt, Laser Communication Systems . New York:. 
John Wiley and Sons, 1968, p. 123.) 




SPECTRAL IRRADIANCE H(X) (Wm 2 M m _ 
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Figure 4.5. Spectral irradiance of brightest stars outside the terrestrial 
atmosphere. (From W. K. Pratt, Laser Communication Systems . 
New York: John Wiley Sons, 1968, p. 126. ) 
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<N>g = expected number of background counts/time slot 

and 

M = number of time slots/word. 

4.4 Calculation of Signal and Background Counts 

In Sections 4.2 and 4.3 we derived the mathematical relationships 
governing the amount of signal and noise counts emitted by the 
photodetector given the signal (laser) and background power present. 
Now we use some typical component parameters to calculate numerical 
results. 

The following parameter values were chosen for a typical 
satellite-to-satellite link: 

transmission wavelength 
transmitter efficiency 
receiver efficiency 
receiver planar field of view 
receiver aperture diameter 
transmitter aperture diameter 
laser power 
quantum efficiency 
data rate 
separation length 
interference filter bandwidth 
Using the above parameter values in (4.11) yields 

<N> max = 77 counts/signal pulse. 


X = 830 nm 
hj, = 0.7 

Uj. * 0.6 

0^=1 mrad 

d = 0.25 m 
r 

tip = 0.1 m 

P L = 200 mW 
n d = 10% 

Q 

D = 2 X 10 words/sec 
z = 3.57 X 10 7 m 
X i = 20 A. 
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This is the maximum number of signal counts we can expect under ideal 
(no pointing error) conditions. Note that pointing error tolerances are 
severe since the Airy intensity pattern falls away rapidly from the 
central maximum. In fact, a pointing error of only 5 prad causes <N> to 
fall to 29 counts/signal pulse. 

Using (4.16) and (4.20) to calculate background counts due to solar 
radiation reflected from the earth filling the receiver field of view 
gives 

<N> ' = 3.08/M counts/time slot (no clouds). 

D 

As previously discussed, if extensive cloud cover were present on earth, 
this value would increase to approximately 30.8/M counts/time slot. 

The moon is the other major source of background radiation. 
Figure 4.4 shows that the spectral irradiance of a full moon is roughly 
10 ^ W/m^/pm. Substituting this value into (4.19) and (4.20), the 
maximum number of background counts possible from a full moon is 

<N> B = 1.52/M counts/time slot. 

This number of counts Is attained only if the receiver is pointing 
directly at a fully illuminated moon. The brightest star with a 
wavelength larger than 0.5 pm has a spectral irradiance of about 
5 X 10 ® W/m^/pm. Inserting this value into (4.17) and (4.20) results 
in 

<N>„ = 6.16 X 10 /M counts/time slot/star , 

B 
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a value small enough to neglect even from many stars of the same 
spectral irradiance. 

Background counts can also occur by dark counts emitted from the 
detector. However, dark counts are usually on the order of 10 to 20 
counts/sec. With the high data rates under consideration, the dark 
counts per time slot are extremely small and also negligible. 

The above calculations show that the worst possible communication 
situation is encountered by a satellite in geosynchronous orbit pointed 
directly towards a cloud-covered section of the earth that is illumi- 
nated by the sun. Using 1 to 3 background counts/time slot for the 
nominal case and 7 to 10 background counts/ time slot for the worst case 
would effectively represent the range of expected operating conditions. 
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CHAPTER 5 

SIMULATION OF THE RECEIVER 


The computer simulation of the receiver is comprised of two main 
parts. The first is the modeling of the photodetector which generates 
the arrival time records, and the second is the simulation of the PLL. 

To generate the arrival time records we follow the procedure 
described in [11]. We know that the interarrival times of a homogeneous 
Poisson process with intensity, X^, are independent and identically 
exponentially distributed also with parameter X^. Using this property 
we first generate a sequence of independent and uniformly distributed 
random variables u^ , , •••, u^ on the interval (0,1). Then using the 
transformation 


== -ln(u jL ) 


(5.1) 


a sequence of independent, exponentially distributed random variables of 
unit parameter is obtained. To generate a homogeneous Poisson process 
we assign l as the arrival time of the first photon, and then define 
as the interarrival time of the second photon and as the interarrival 
time of the photon. Now, the arrival time of the k C ^ photon, x^, is 



(5.2) 
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Next a random PPM intensity sequence, ^(t), is created. The final 
arrival time sequence used in (2.5) is formed by solving 



(5.3) 


for the t^’s [11]. Appendix A.l contains a listing of the Fortran pro- 
gram ARRIVAL used to generate the arrival time record for an M-ary PPM 
signal using the described method. 

The PLL was simulated using the Advanced Continuous Simulation 
Language (ACSL). This language is designed for modeling and evaluating 
the performance of continuous systems that are described by time 
dependent, nonlinear differential equations and transfer functions. 
ACSL allows direct transformation of the mathematical diagram of 
Figure 3.2 into a suitable simulation model program. A listing of the 
ACSL simulation program for the PLL model is in Appendix A. 2. This 
program takes an externally prepared arrival time sequence, processes it 
through the filter and squaring loop, constructs the shot noise, applies 
the resulting signal to the PLL, and gathers data on the phase error. 
The simplicity of ACSL is evident in lines 30-35 of the simulation 
program. These six lines contain all the model information necessary to 
implement the PLL. In fact, the filter is completely described on line 
33 by the ACSL statement REALPL which implements a simple low-pass 
filter. The integrating effect of the VCO is contained in line 35 by 
use of the INTEG statement which performs an integration on the VCO 
input control voltage. Most of the remaining program is devoted to 
simulation of the preprocessing components and shot noise construction. 
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To show the validity of the simulation model, we present results of 
two simulations: one without signal preprocessing and the second with 
preprocessing. From the results in Chapter 2, we expect phase-lock 
behavior only for the preprocessed case due to the presence of the 
dominating spectral frequency component in the PLL input. Figures 5.1 
and 5.2 show the Fourier transforms of the PLL input and VCO output 
signal when preprocessing is by-passed, and where the frequency scale is 
normalized to the time slot rate. Note the lack of a dominating 
frequency in the input spectrum and that the VCO output spectral 
component occurs below the slot rate at , the VCO free-running 
frequency. Figure 5.3 is a plot of Ao)t versus time. Recall that 
Aoit=K o v c ( t)t is the phase change of the VCO due to the control voltage 
applied to its input. This change, when added with the phase due to the 
free-running VCO frequency, should equal the phase of the input signal. 
The frequency difference between the slot rate and was 0.103 rad/sec, 
which should equal the slope of Figure 5.3. Instead, Figure 5.3 shows a 
very random behavior rather than a well-defined slope of 0.1 rad/sec. 
Clearly with no signal preprocessing, and thus no dominating spectral 
component, the loop cannot achieve phase lock. 

Figures 5.4 and 5.5 display the Fourier transforms of the input and 
VCO output when preprocessing is used. Again, the frequency scales are 
normalized to the slot rate. Notice the very strong slot frequency 
component in the PLL input signal. Now the VCO output spectrum in 
Figure 5.5 has shifted from its free-running frequency to match the PLL 
input signal . frequency component at the slot rate. The plot of Acut 
versus time of Figure 5.6 now has the desired slope of 0.1 rad/sec. 
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NORMALIZED FREQUENCY (cu/w s ) 


Figure 5.1. Phase-locked loop input signal frequency spectrum when no 
preprocessing is used. Note the lack of a dominating fre- 
quency component. 
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Figure 5.2. VCO output frequency spectrum when no signal preprocessing is 
used. Note that the dominating frequency component occurs 
below the slot rate at the VCO free running frequency. 
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NORMALIZED FREQUENCY (cl»/oj s ) 


Figure 5.4. Phase-locked loop input signal frequency spectrum when 

preprocessing is used. The dominating spectral component 
occurs at the slot frequency. 
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Figure 5.5. VCO output spectrum when preprocessing is used. Note that 
the spectral component has shifted up to the slot frequency 
from its free running frequency. 
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Figure 5.6. Aoot versus time when preprocessing is used. Note the slope of 
approximately 0.1 rad/sec indicates tracking performance of 
the phase-locked loop. 
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CHAPTER 6 

COMPUTER SIMULATION RESULTS 

Computer simulation of the timing subsystem was conducted under 
varying signal and background strengths to represent a wide range of 
operating conditions. In this chapter we present the simulation results 
and compare them to a current theoretical expression. 

6.1 Simulation Performance Parameters 

The parameter of interest is the phase error, <j>(t), defined as 

♦<t) = e s (t) - 8 o (t) (6.1) 

where, again, 0 g (t) and 0 q C t ) are the PLL input and VCO output signal 

phases, respectively. Once the loop is in lock, we expect the phase 

error to vary about some constant value. Assuming <j>(t) is approximately 

Gaussian, the best system performance indicator of the PLL is the RMS 

phase error, a, (rad). 

9 

; Derivation of an exact expression for a, is lengthy and difficult. 

9 

However, C. Chen has formulated the following approximation of o for 

9 

the receiver timing subsystem presented above. From [12], 
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where 

K g = expected signal count 

H l (<d) * closed-loop frequency response of the PLL 

P(u>) = Fourier transform of received pulse shape 

H( w) = frequency response of the detector and linear filter 

M = number of time slots/word 

T = slot width, 

s 

An important feature of (6.2) is that a ^ is inversely proportional to 

the square root of the average signal count, K . This important 

s 

property provides a convenient means of comparing the theory to the 
actual simulation behavior in the next section. 


6.2 Numerical Results 

A total of 26 computer simulations were completed under varying 
signal and background count conditions for both 4 slot and 8 slot PPM 
signals. The following is a summary of system parameters used in the 
trials: 


time slot width 
time slot rate 
free running frequency 
phase detector gain 


T =1 sec 
s 

a) = 2 ir rad/ sec 
s 

Wf = 6.18 rad/sec 

K, = 0.3333 V/rad 
a 
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VCO gain constant K q = 1 rad/sec/V 

loop filter bandwidth 1/t = 0.3333 rad/sec. 

Each trial represents an 8000 sec operation time, corresponding to 2000 

data words (4000 bits) for 4 slot PPM and 1000 data words (3000 bits) 

for 8 slot PPM. The phase error was calculated at the beginning of each 

time slot, and these values were used to calculate the RMS phase error 

over the simulation time. Plots of a, versus time for each trial appear 

9 

in Appendix B. Notice as the average background count, K^, increases 

and the average signal count, K g , decreases, the RMS phase error 

increases and exhibits progressively erratic behavior. 

The versus time plots of Appendix B are summarized in Table 6.1. 

Table 6.1 lists K^, K g , o^, as well as the average signal count/bit, 

Kfait , so 4 and 8 3lot PPM formats can be correctly compared at the bit 

level. On a signal count/bit comparision, 8 slot PPM has a higher a, 

9 

(and thus inferior performance) than 4 slot PPM with the same and 

K fa . Figure 6.1 shows that the difference in performance increases as 
Kfait increases. Also, two of the simulation trials for 8 slot PPM 
(K g =45, K^=5 , and K s =90, K^=20) exhibited cycle-slipping behavior, 

indicating that 8 slot PPM may be more sensitive to input signal 
variations. 


Figure 6.2 is a log-log plot of a, versus average signal count for 4 

9 

slot PPM and a constant of 1, 5, and 10, respectively. Figure 6.3 is 

a similar plot for 8 slot PPM with K^=5 only. As indicated in 

Section 6.1, we expect o to be inversely proportional to the square 

9 

root of K g , corresponding to an ideal slope of -0.5 in Figures 6.2 and 
6.3. Each of the lines drawn in Figures 6.2 and 6.3 represents a least 
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NUMERICAL RESULTS OF SIMULATION TRIALS 


Average 

Signal 

Count 

K 

s 

Average 

Signal 

Count/Bit 

K bit 

PPM 

Signaling 

Order 

M 

Average 

Background 

Count 

*b 

RMS 

Phase 

Error 

o, (rad) 
<P 

30 

15 

4 

1 

0.2926 

30 

15 

4 

5 

0.3399 

40 

20 

4 

1 

0.2609 

40 

20 

4 

5 

0.2737 

40 

20 

4 

10 

0.3146 

50 

25 

4 

1 

0.2496 

50 

25 

4 

5 

0.2303 

50 

25 

4 

10 

0.2572 

60 

30 

4 

0 

0.2181 

60 

30 

4 

1 

0.2235 

60 

30 

4 

5 

0.2239 

60 

30 

4 

10 

0.2402 

60 

30 

4 

20 

0.2594 

80 

40 

4 

1 

0.1974 

80 

40 

4 

5 

0.1857 

80 

40 

4 

10 

0.2070 

45 

15 

8 

5 

0.4257 

60 

20 

8 

5 

0.3573 

75 

25 

8 

5 

0.3234 

80 

26.67 

8 

5 

0.3217 

90 

30 

8 

0 

0.2990 

90 

30 

8 

1 

0.2978 

90 

30 

8 

5 

0.3029 

90 

30 

8 

10 

0.3053 

90 

30 

8 

20 

0.3823 

120 

40 

8 

5 

0.2805 
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squares fit of data points with respect to a constant value of K, . In 

D 

Figure 6.2 the slopes for 4 slot PPM trials are: -0.40 for K,=l, -0.60 

for 1^=5, and -0.58 for K^IO. For 8 slot PPM trials in Figure 6.3, the 

slope was -0.425 for 1^=5. These results compare favorably to the ideal 

slope of -0.5 predicted by (6.2). Currently, a factor of between three 

and four separates the theoretical and simulation values. It is 

reemphasized that (6.2) is an approximation, not an exact expresssion. 

The . exact expression for a, is unwieldy and difficult to evaluate. 

9 

Efforts are continuing to obtain a more accurate approximation of the 

exact expression for a, [12]. 

9 

Figure 6.4 shows the relationship, for 4 slot PPM, between a. and K, 

(J> D 

when K is held constant to 40, 50, 60, and 80 average signal counts, 
s 

Figure 6.5 is a similar plot for 8 slot PPM for K =90 only. In all 

. s 

cases, the relationship between the RMS phase error and background count 

is not a simple one. For lower background counts, the RMS phase error 

rises moderately with increasing K^. As the background count continues 

to increase, the RMS phase error rises sharply at about 1^=5, Indicating 

a rapid deterioration of timing performance. 

The ratio between the timing offset, AT, and the time slot width, 

T , called the time slot error (TSE), is related to the RMS phase error 
s 


XSE = — = £ 

TS T 2 TT 
s 


(6.3) 


An increasing TSE indicates receiver timing degradation and an 
increasing probability of bit error, PBE. In fact, the advantages 
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gained in block coding (when the probability of word error (PWE) is 
converted to PBE at normalized rates) are almost eliminated as the 
timing error (or TSE) increases [4j. As a typical example, consider the 
4 slot simulation run of K g =60 and K^=l of Table 6.1. The RMS phase 
error is 0.2235 rad/sec, corresponding to a TSE of 0.0356. This 
indicates the receiver timing is off by 3.56%. The 8 slot simulation 
trial of K g =90 and K^=l has a RMS phase error of 0.2978 rad/sec and a 
TSE of 0.0474, corresponding to a 4.74% timing error. Detailed 
discussion and analysis of timing error effects on optical PPM 
communication systems can be found in [1] and [4]. 

This concludes the numerical results of the computer simulation 


data 
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CHAPTER 7 
CONCLUSION 


In this report we have investigated the timing performance of a 
phase-locked loop (PLL) in an optical pulse position modulation (PPM) 
communication system. It was shown that the output from a photodetector 
alone contains no dominating spectral frequency component and is 
insufficient as a signal source for the PLL. A simple filter and 
squaring loop was shown to generate a strong spectral component at the 
slot frequency. Preprocessing the signal in this manner enables the PLL 
to lock up, track the signal, and provide timing information to the 
receiver decision section. 

We have discussed various PLL design considerations. Once design 
parameters are fixed, a change in the input signal magnitude can alter 
the PLL performance. To alleviate this problem, an expression for an 
attenuation constant was derived which would limit the average PLL input 
signal to unity, allowing the intended design parameters to remain 
unchanged over a wide range of signal inputs. 

A general satellite-to-satellite link equation was derived to 
calculate the average number of signal counts emitted by the receiver 
detector. Average signal count depends on many factors, including 
physical component values (aperture size, wavelength, etc.), data rate, 
and pointing error. Sources of background radiation such as the moon. 
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stars, and earth were examined and it was shown that a cloudy, 
illuminated section of the earth results in the largest amount of 
background radiation at the receiver. Typical component values for the 
link were chosen, and average signal and background counts calculated. 

Finally, the entire timing subsystem was simulated on a digital 
computer and the JIMS phase error of the PLL was calculated. The 
simulations showed the RMS phase error is (approximately) inversely 
proportional to the square root of the average signal count for both 4 
slot and 8 slot PPM signals. This agrees with a theoretical expression 
recently derived for this timing subsystem which predicts an exact 
inverse relationship between the RMS phase error and the square root of 
the average signal count. 
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APPENDIX A 
PROGRAM LISTINGS 


A.l Arrival 


1 

2 C 

3 C 

4 C 

5 C 

6 C 

7 C 

8 C 

9 C 

10 C 

11 C 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 10 
22 

23 

24 100 

25 

26 170 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 180 

37 200 

38 

39 

40 

41 

42 C 

43 C 


PROGRAM ARRIVAL 

THIS PROGRAM CONSTRUCTS THE PHOTON ARRIVAL 
TIMES FOR A M-ARY PPM DATA SEQUENCE DETECTED 
BY A PKOTODETECTOR . 

PULSHI=AMPLITUDE OF SIGNAL PULSE 
NSLOT=ORDER OF THE PPM SIGNAL 
TMAX=LENGTH OF RECORD DESIRED IN SECONDS 
DARR=BACKROUND RADIATION PER PULSE 
TSLOT=DURATION OF A SINGLE PULSE 

DIMENSION ARIVAL( 1024) 

INTEGER ISLOT( 25) 

DATA NSLOT/ 4/, TMAX/ 9000.0/ 

DATA PULSHI/40 . / , DARK/ 10./, TSLOT/ 1 .0/ 

AREA =0.0 
TNOW = 0.0 
TIMLST =0.0 
CALL RANSET(31) 

DO 10 1=1,10 

U = RANFC) 

CALL DATGEN(ISLOT, NSLOT) 

KK = 1 

CALL POISSN(ARIVAL, TIMLST, 1024, 1.0) 

DO 200 I = 1,1024 

DIFF1 = (ISLOT(KK) * PULSHI + DARK) * TSLOT 
DIFF2 = ARIVAL(I) - AREA 
IF (DIFF2 .LE. DIFF1) GO TO 180 
KK = KK+1 

TNOW = TNOW + TSLOT 
AREA = AREA + DIFF1 
IF (KK.LE. NSLOT) GO TO 170 
CALL DATGEN( I SLOT, NSLOT) 

KK = 1 
GO TO 170 

ARIVAL(I) = TNOW + DIFF2/DIFF1 

CONTINUE 

WRITE(3) (ARIVAL(IW) , IW=1 , 1024) 

IF (TNOW .LT. TMAX) GO TO 100 

RETURN 

END 
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C 

45 


46 

C 

47 

C 

48 

C 

49 

C 

50 

C 

51 


52 

1 

53 


54 


55 


56 

12 

57 


58 


59 

10 

60 


61 


62 


63 

C 

64 

C 

65 


66 

c 

67 

c 

68 

c 

69 

c 

70 


71 


72 


73 


74 

10 

75 


76 


77 


78 


79 



SUBROUTINE POISSN( ARIVAL , TIMLST , NUMTRY , PARAM) 

SUBROUTINE POISSN GENERATES THE POISSON INPUT 
SEQUENCE OF UNIT PARAMETER USING UNIFORM 
RANDOM VARIABLES. 

DIMENSION ARIVAL(l) 

UNIF = RANF() 

IF (UNIF. EQ. 0.0) GO TO 1 

ARIVAL ( 1 ) = TIMLST - <AL0G( UNIF) /PARAM) 

DO 10 J=2, NUMTRY 

UNIF = RANFQ 

IF (UNIF. EQ. 0.0) GO TO 12 

ARIVAL(J) = ARIVAL (J-l ) - (ALOG( UNIF) /PARAM) 

CONTINUE 

TIMLST = ARIVAL( NUMTRY) 

RETURN 

END 


SUBROUTINE DATGEN ( I SLOT , NSLOT ) 

SUBROUTINE DATGEN GENERATES THE RANDOM 
SIMULATED PPM DATA SEQUENCE. 

INTEGER ISLOT(l) 

INTEGER MASK, IRAND 
MASK = NSLOT 
DO 10 1=1, NSLOT 

ISLOT(I) = 0 

IRAND = IFIX(RANF( ) * 1000) 

IRAND = MOD( IRAND, MASK) 

I SLOT ( IRAND+1 ) = 1 

RETURN 

END 
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A. 2 Simulation of a Phased-Locked Loop 


1 PROGRAM SIMULATION OF A PHASE-LOCKED LOOP 

2 "THIS PROGRAM IS DONE IN THE ADVANCED" 

3 "CONTINUOUS SIMULATION LANGUAGE (ACSL)" 

4 •• •• 

5 INTEGER INDEX 

6 REAL W0,KV 

7 REAL ARIVAL(2049) 

8 CONSTANT THETA = 0..KV =.3333,... 

9 FIC2 = 0.0, 

10 PHIO = 0.0, WO =6.18, ... 

11 INITIAL 

12 CINTERVAL CINT = 0.125 

13 NSTEPS NSTEP = 3 

14 IALG =4 

15 INDEX =1 

16 ENDTIM = 8050.0 

17 PROCEDURAL ( ARIVAL=ARIVAL) 

18 CALL SETVAL(ARIVAL) 

19 END $"OF PROCEDURAL" 

20 END $ " OF INITIAL " 

21 " " 

22 DYNAMIC 

23 DERIVATIVE $"THIS SECTION CONTAINS THE PLL" 

24 "SIMULATION MODEL EQUATIONS." 

25 • 

26 

27 PROCEDURAL (SHOT=T,ARIVAL, INDEX) 

28 CALL SHOTNZ(T,ARIVAL, INDEX, SHOT) 

29 END$ " OF PROCEDURE " 

30 PHASE = WO*T+PHI+THETA 

31 COPHAS=COS( PHASE) 

32 MIXOUT = SHOT*COS( PHASE) 

33 FILOUT = REALPL( 3. .MIXOUT, 0.0) 

34 PHIDOT = KV* FILOUT 

35 PHI = INTEG( PHIDOT, PHIO) 

36 TERMT(T.GE. ENDTIM) 

37 END $ " OF DERIVATIVE " 

38 WRITE(l) PHASE $ "SAVES THE VCO OUTPUT PHASE FOR ANALYSIS" 

39 END $ " OF DYNAMIC SECTION " 

40 END $ " OF PROGRAM " 

41 

42 



t 


71 


43 

44 C 

45 C 

46 C 

47 C 

48 C 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 C 

59 C 

60 C 

61 C 

62 C 

63 C 

64 

65 

66 

67 10 

68 

69 

70 

71 20 

72 

73 

74 

75 100 

76 

77 

78 101 

79 

80 

81 C 

82 C 

83 C 

84 C 

85 C 


SUBROUTINE SETVAL(ARIVAL) 

SETVAL READS IN THE FIRST TWO BLOCKS (2048 VALUES) OF 
THE ARRIVAL TIMES FROM THE EXTERNALLY PREPARED ARRIVAL 
TAPE. 

REAL ARIVAL(2049) 

READ(3) (ARIVAL( J) , J=1 , 1024) 

R£AD(3) (ARIVAL(J) ,J=1025,2048) 

RETURN 

END 


SUBROUTINE SHOTNZ (T.ARIVAL, INDEX, SHOT) 

SUBROUTINE SHOTNZ EVALUATES THE SHOTNZ SIGNAL 
FOR INPUT TO THE PLL AT A GIVEN TIME T. SHOTNZ 
RETURNS THE CURRENT VALUE IN THE PARAMETER SHOT. 


REAL ARIVAL(2048) 

INTEGER INDEX 
Q=T+10 

IF (ARIVAL( INDEX ).GE.Q) GO TO 100 
INDEX = INDEX +1 
IF (INDEX.LT. 2049) GO TO 10 
DO 20 1*1,1024 

ARIVAL(I) * ARIVAL( 1+1024) 

READ(3) (ARIVAL(I), 1=1025, 2048) 

INDEX = 1025 
GO TO 10 
SHOT =0.0 
FLAG =1 
II = 0 
II =11+1 

SHOT = SHOT + GATE( Q-ARIVAL( INDEX-II) , IFLAG) 
IF ( IFLAG. GE.O) GO TO 101 

THE NEXT LINE TAKES THE SHOT NOISE VALUE AND 
SQUARES IT AND DIVIDES BY THE ATTENUATION 
CONSTANT OF THE PREPROCESSING SECTION. 
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87 

88 

89 

90 C 

91 C 

92 

93 C 

94 C 

95 C 

96 C 

97 

98 

99 
100 
101 
102 

103 

104 10 

105 

106 
107 


SH0T1 = SHOT*SHOT/40. 528 
SHOT = SH0T1 

RETURN 

END 


REAL FUNCTION GATE (X, IFLAG) 

FUNCTION GATE SIMULATES THE COMBINED 
PREPROCESSING FILTER IMPULSE RESPONSE. 

REAL PI2 

DATA PI2/6. 283185307/ 

IF (X.GT.1.0) GO TO 10 
G = SIN(PI2*X) 

IFLAG « 1 
GATE = G 
RETURN 

IFLAG ® -1 
GATE - 0.0 
RETURN 
END 



RMS PHASE ERROR (RAD) 
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APPENDIX B 

RMS PHASE ERROR VERSUS TIME FOR SIMULATION TRIALS 



~~ SECONDS 

Figure B.l. Simulation trial for K = 30, K, = 

s b 


1 and M = 4 



RMS PHASE ERROR (RAD) 
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SECONDS 


Figure B.2. Simulation trial for K g = 30, = 5 and M = 4. 


RMS PHASE ERROR (RAD) 
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SECONDS 

Figure B.3. Simulation trial for K g = 40, = 1 and M = 4. 



RMS PHASE ERROR (RAD) 


76 



1000. 5000. 5000. 7000. 

SECONDS 

Figure B.4. Simulation trial for K = 40, K, = 5 and M = 4. 

s b 


RMS PHASE ERROR (RAD) 



1000 . 5000 . 5000 . 7000 . 

r SECONDS 
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1000. 5030. S303 . 7333. 

SECONDS 

Figure B.6. Simulation trial for K g = 50, = 1 and M = 4. 
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1000. 5000. 5000. 7000. 

SECONDS 

Figure B.7. Simulation trial for K = 50, K, = 5 and M = 4. 

s b 


RMS PHASE ERROR (RAD) 
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SECONDS 

Figure B.8. Simulation trial for K = 50, K, = 10 and M = 4. 

s D 


RMS PHASE ERROR (RAD) 
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1030. 5000. 5000. 7000. 

V • SECONDS 

Figure B.9. Simulation trial for K g = 60, = 0 and M = 4. 
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Figure B.13. 


5000. 5000. 7000 

SECONDS 


Simulation trial for K = 60, K. = 20 and H 

ft 9 n 
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1000. 5000. 5000. 7000. 

SECONDS 

Figure B.19. Simulation trial for K g = 75, = 5 and M = 8. 








2000 : 
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Figure B.25. 


5000. 5000. 7000. 

SECONDS 


Simulation trial for K 


90, = 20 and M = 8. 
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